Accelerating an adiabatic process by nonlinear sweeping 
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We investigate the acceleration of an adiabatic process with the same survival probability of the 
ground state by sweeping a parameter nonlinearly, fast in the wide gap region and slow in the 
narrow gap region, as contrast to the usual linear sweeping. We find the expected acceleration in 
the Laudau-Zener tunneling model and in the adiabatic quantum computing model for factorizing 
the number N = 21. 
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Due to its robustness against noises, adiabatic quan- 
tum computing (AQC) is one of the promising quan- 
tum computation schemes which are experimentally fea- 
sible^. However, the main concern is that the energy 
gap between the ground state and the lowest excited state 
of a quantum computer decreases exponentially with the 
increasing number of qubits 3 '-, so that the adiabatic evo- 
lution time required by the algorithm would increase ex- 
ponentially and eventually cancels the speedup benefit of 
a quantum computer. 

One way to overcome the long-adiabatic-evolution- 
time obstacle is using rapid adiabatic passage method^— , 
which is capable of transferring an initial ground state to 
a final ground state within an arbitrary short time. The 
price is that all the gap information must be provided, 
which is usually beyond the reach of current techniques 
for a many-qubit quantum computer. It is thus useful to 
accelerate an adiabatic process, similar to the rapid adi- 
abatic passage method, but with less information about 
the gap, for instance, only the extremes of the gap during 
the total evolution. 

In this paper, we investigate nonlinear sweeping meth- 
ods to accelerate adiabatic processes, including the fa- 
mous Laudau-Zener (LZ) tunneling model and the AQC 
model 9 '-—, with only the knowledge of the extremes of 
the gap. For the LZ model, we numerically confirm that 
the acceleration effect of a nonlinear sweeping. We then 
develop a similar nonlinear sweeping scheme for a simpli- 
fied AQC model and finally confirm the effectiveness of 
this nonlinear sweeping method by factorizing the num- 
ber N = 21. 

The LZ model of a two-level system is described by 



H(t) = huj x o~x + hio z {t)a z 



(1) 



where uj XyZ is the frequency of the Zeeman splitting along 
the x or z direction, and a x<z the Pauli matrix. For a 
symmetric system as shown in the left panel of Fig. [TJ 
the time dependence of a linear sweeping of co z (t) is 



w,(t) 



-U)Q + Vt, 



(2) 



where v = 2loq/T denotes the sweeping velocity, with 
T being the total evolution time. If w z (t) sweeps fast 
from —ujq to wq, an initial ground state may have a large 
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FIG. 1: The instantaneous energy level diagram for a two- 
level LZ model (left) and the fidelity of the final ground state 
(right). The parameters are cjo = 10, Wi = 1, and T = 10 
for the left panel. We set h = 1 for convenience. In the 
right panel, the dash-dotted line denotes the fidelity under the 
linear sweeping and the dashed line the nonlinear sweeping. 
The fidelity for the nonlinear sweeping is higher than that for 
the linear sweeping. 



probability to be transferred to the final excited state. 
Such a probability is described by the Laudau-Zener for- 
mula^. If ui z (t) slowly changes, the system would stay in 
the ground state. 

It is for sure that the linear sweeping of ui z {t) is not 
optimized for the fidelity F , which measures the survival 
probability of the final ground state and is defined as 



F=|WT)|V 3 )| i 



(3) 



with \tp g ) denoting the instantaneous ground state of the 
Hamiltonian H = H(T) and \ip{T)) the actually evolved 
wave function at the end time t = T. An intuition for the 
optimization is sweeping fast if the energy gap is large but 
sweeping slow if the gap is small, where the energy gap 
A is defined as the energy difference between the ground 
state uig and the lowest excited state w e , i.e., A = uj e ~u g . 
It is straightforward to construct a nonlinear sweeping 
form which satisfies the above intuition, 



Wz(t) 



-OJ 



to 



t<t c 
t > t. 



(4) 
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with t c = T/2. The gap A is minimal at t c 
sponding sweep velocity is 



v(t) = 2wq 



t 



tc 



- 1 



The corre- 



(5) 



where 



which is large near t — and T where the gap is large but 
small near t = t c where the gap is small. The velocity 
for the whole evolution can be adjusted by changing the 
total time T. The longer the T is, the smaller the velocity 
is. 

We present in the right panel of Fig. [T] the fidelity F 
of the hnal ground state under either the linear sweeping 
form Eq. @ or the nonlinear sweeping form Eq. (Q}. It 
clearly shows that the fidelity with nonlinear sweeping 
is much higher than that with linear sweeping, provided 
that the total evolution time T is the same. In other 
words, the adiabatic LZ process can be accelerated by 
adopting a nonlinear sweeping form. 

We note that other more complex nonlinear forms also 
accelerate the adiabatic process, especially the stimu- 
lated Raman adiabatic passage or the hyperbolic tangent 
form££, and the fidelity F may be even higher. These 
methods usually require all the information of the en- 
ergy levels which is not always available. While in the 
nonlinear formula, only the position of the minimal gap t c 
is required. In this sense, the nonlinear sweeping method 
is easier to implement in experiments. 

We consider next a time-dependent two-level system, 
that is a single qubit AQC model, 



H(t) = [1 - s(t)]uj x a x + s(t)uj z a z , 



(6) 



where s(t) is increasing monotonically with s(0) = and 
s(T) = 1. By sweeping slowly s(t) from to 1, the two- 
level system transfers from an initial ground state to its 
final ground state. A linear sweeping and a quadratic 
sweeping form used in experimenti^ are, respectively, 



8{t) 



t 

f 



and 



*(*) = ( t 



(7) 



(8) 



It is obvious that neither the linear sweeping nor the 
quadratic sweeping satisfies the requirement of fast 
sweeping in large gap region and slow sweeping in small 
gap region. We find the optimal sweeping form as follows. 

It is the best to do the similar procedure to the above 
LZ model and utilize the optimized sweeping scheme. We 
alternate the single qubit AQC Hamiltonian Eq. (|6|) to a 
LZ model Eq. (JTJ) where only one of the two parameters 
before the Pauli matrices is time-dependent 



H(t) = ujj_a^ + uj n (s(t))a n , 



(9) 



a± = a z sm v + a x cos ( 

o"n = o z cos 6 — o x sin ( 

cjj_ = uj x cos 9, 

Unit) = s(t)fl — OJxWO.6 



with cos0 = w z /f2, sva.9 = uj x /Q,, and SI = \/ui% -f 
To locate the position of the minimal gap at s(t c 
we calculate the energy gap 



A(a) = + (1 



where the minimum satisfies 

dA 
ds 

It is easy to find that 



= 0. 



UJ, 



(10) 



(11) 



(12) 



Another way to obtain the above result of s c is setting 
w n (s c ) = 0. 

To utilize the optimized quadratic sweeping scheme, 
for the LZ model we notice that s in the AQC model 
Eq. (J9j) takes the place of t in the LZ model Eq. (QJ, so 
we use the following relation 



duj n (s(t)) 
dt 



A\s(t) 



(13) 



The value of A can be determined by the boundary of 
w n (s(t)). This equation implies that 



ds 
~dl 



(14) 



which gives the relation s(t). Similar to the optimized 
scheme in the LZ model, the form of s(t) sweeps fast in 
the large gap region and slow in the small gap region with 
respect to s. 

It is straightforward to obtain the relation s(t), which 
satisfies not only the initial and final values but also an 
extreme point s(t c ) = s c , 



s(t) = 



-(1 



-at/t c ^ 



1 



a OLt/t c 



t < t c 
t > t c 



(15) 



We have introduced an adjustable parameter a varies 
the curvature of s(t) [see the left panel of Fig. [2]. When 
a = 0, s(t) is linear; When a = 1, s(t) is a two-interval 
exponential-like piecewise function; When a — ¥ oo, s(t) is 
step-like function. To determine t c , we adopt a simplest 
form t c /T = s c . 

We present in the right panel of Fig. [5] the fidelity of 
the final ground state under the time-dependent Hamil- 
tonian Eq. ^) for various sweeping schemes. Note that 
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FIG. 2: Left: Various forms of s(t); Right: Dependence of 
the fidelity on the total evolution time. Dotted line with 
circles — Exponential-like with opitimized a; Dotted lines — 
Exponential-like with a — 5; Solid lines — Exponential-like 
with a — 1; Dash-dotted lines — Linear; Dashed lines - 
Quadratic. The parameters are ui x = 18 and u) z — 30. The 
Exponential-like schemes are better than the Linear of the 
Quadratic one. 




S T 

FIG. 3: Left: The instantaneous energy gap versus s for 
N — 21. The minimal gap position is near s c ~ 0.74. 
Right: The dependence of the fidelity on the total evolution 
time for factorizing N with optimized exponential-like sweep- 
ing scheme (dotted line with circles), linear sweeping scheme 
(dash-dotted line), and quadratic sweeping scheme (dashed 
line). The parameter is g = 30. Among the three schemes, 
the optimized exponential-like one shows the highest fidelity 
at the same evolution time T and is the best choice. 



the initial state of the two-level system is an eigenstate 
of a x , which overlaps with the final instantaneous ground 
state so the fidelity is 0.5 even if T approaches zero. As 
shown in the right panel of Fig. [3J it is obvious that the 
quadratic scheme is the worst choice. The fidelity with 
the linear scheme is between that with quadratic scheme 
and that with exponential-like scheme. By optimizing a 
for each T in the exponential-like schemes, the fidelity 
reaches its highest value. We also notice that for some 
values of a in the exponential-like schemes, the fidelity 
oscillates if the total evolution time T becomes large, 
e.g., a = 5. These oscillations might be due to the over- 
long-time stay of the system near the minimal gap region 
where the gap is almost constant^. 

The time-dependent Hamiltonian for factorizing N = 
21 is an interpolation of an easy-realized initial Hamilto- 
nian Hq and a problem-solved Hamiltonian Hp^, 

H(t) = [1 - s(t)] H + s(t) Hp, (16) 

where s(t) changes from to 1 as the system evolves and 

Hq = g{(Jlx + &2x + <?3x), 

H P = [N-(2I-a lz ){AI-a 2z -2a 3z )} 2 . 

In the above equations, g is the Zeeman splitting of a 
uniform magnetic field along x direction, I is the 2-by-2 
identity matrix, and <7j a with i — 1,2,3 and a — x, z 
denotes the a Pauli matrix for the ith spin. The ground 
state of Hp, | 4-i) ® | 4-24-3) interpreted as 3 ® 7, is the so- 
lution to the factors of N. The 3-spin system is prepared 
initially in the ground state of Hq. After an adiabatic 
evolution under the total Hamiltonian Hit), the system 
eventually reaches the ground state of H(T) = Hp and 
the number N is factorized. 

Usually, two sweeping schemes are experimentally 
adoptedii^. One is the linear scheme as defined in 



Eq. (J7J and the other is the quadratic scheme as in 
Eq. ((SJ). We have shown in the two-level model that both 
the linear and the quadratic sweeping schemes are not 
the optimal. We may use the exponential-like scheme to 
accelerate the AQC process if the energy gap between 
the ground state and the lowest excited state is approx- 
imately in an inverted parabolic shape with respect to 
s. Fortunately, the gap for N = 21 is indeed very sim- 
ilar to that for the single qubit AQC model system, as 
shown in the left panel of Fig. [3J Thus we directly ap- 
ply the exponential-like sweeping scheme to the adiabatic 
process N — 21 factorization. 

We present in the right panel of Fig. [3] the fidelities 
of the evolved final state to the ground state of Hp at 
different total evolution time T with the linear sweeping 
scheme, the quadratic scheme, and the exponential-like 
scheme. We limit ourselves in the short T region where 
the difference among these schemes are prominent. It 
is clearly shown in the figure that the fidelity with the 
exponential-like sweeping scheme is higher than that with 
either the linear or the quadratic sweeping scheme for the 
same total evolution time T. Thus, to reach the same 
final fidelity, the exponential-like sweeping scheme uses 
the shortest time and accelerates the AQC process. 

For larger N, more qubits are required to factorize the 
number. The energy levels inevitably become so complex 
that the energy gap A may show multiple extremes. In 
this situation, a more complecated piecewise function s(t) 
with exponential-like form in each interval might be more 
appropriate. More explorations along this direction is 
worthy in the future. 

In conclusion, we develop a nonlinear sweeping scheme 
for an adiabatic process, which changes fast in the large 
gap region and slow in the small gap region. This scheme 
improves the survival probability of the final ground 
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state and may thus accelerate the adiabatic process. We 
confirm the usefulness of this nonlinear scheme in the 
Laudau-Zener tunneling model, the single qubit AQC 
model, and the AQC of factorizing the number N = 21. 
The proposed nonlinear scheme may be useful in AQC 
experiments to shorten the adiabatic evolution time. 
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